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ABSTRACT 



Context. For relativistic modelling of high-accuracy astronomical data several time scales are used: barycentric and geocentric coordinate 
times, TCB and TCG, as well as two additional time scales, TDB and TT that are defined as linear functions of TCB and TCG, respectively. 
Aims. The paper is devoted to a concise but still detailed explanation of the reasons and the implications of the relativistic scalings of 
astronomical quantities induced by the time scales TDB and TT. 

Methods. We consequently distinguish between quantities and their numerical values expressed in some units. 

Results. It is argued that the scaled time scales, the scaled spatial coordinates and the scaled masses should be considered as distinct quantities 
which themselves can be expressed in any units, and not as numerical values of the same quantities expressed in some difi"erent, non-SI 
units {''TDB units" and 'TT units"). Along the same lines of argumentation the system of astronomical units is discussed in the relativistic 
framework. The whole freedom in the definitions of the systems of astronomical units for TCB and TDB is demonstrated. A number of 
possible ways to freeze the freedom are shown and discussed. It is argued that in the future one should think about converting AU into a 
defined quantity by fixing its value in SI meters. 
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• rH 1. Introduction 

, It is well known that the accuracy of modern astronomical ob- 
■ - - ' servations has attained a level where numerous relativistic ef- 
fects can no longer be ignored. Moreover, the whole set of 
astronomical concepts used for interpretation of observational 
data has to be formulated in the framework of general relativ- 
ity. In recent years significant progress has been achieved in 
this direction. A rigorous post-Newtonian framework for rel- 
ativistic data modelli ng has been adop ted by the International 
Astronomical Union (" Soffel et allifflB . and reference therein). 
Nevertheless, the situation is not yet fully satisfactory. One of 
the main factors retarding the adoption of a fully self-consistent 
relativistic framework for fundamental astronomy is the exis- 
tence of "inertia" or "traditions" which are quite difficult to 
overcome. Some of these traditions are heavily based on spe- 
cial approximations in the framework of Newtonian physics, 
some other are based on a Newtonian-like interpretation of the 
theory of relativity. 

One of the controversial questions of the latter kind is the 
situation with the linear scaling of astronomical time scales and 
spatial coordinates related to the theory of relativity. Although 
this question is clear and even almost trivial from the theoreti- 
cal point of view, practical implications of the scaling are some- 
times tricky and often understood in a confusing way. The aim 



of this paper is to provide a concise, self-consistent and rigor- 
ous description of the whole situation with relativistic scalings. 
Interestingly, the same discussion can be used to clarify the 
definition of the system of astronomical units in the relativistic 
framework. This subject has been only marginally discussed in 
the literature and not all what was published on this subject was 
correct. 

In Section |2] the relations between quantities and their nu- 
merical values expressed in some units are summarized. The 
justification for and implications of TDB, being a scaled ver- 
sion of the coordinate time TCB of the Barycentric Celestial 
Reference System (BCRS), are discussed in Section[3] Section 
|4] is devoted to the relativistic scahng in the Geocentric 
Celestial Reference System (OCRS). The concept of coordi- 
nate time scales is elucidated in Section |5] The difficulties 
that appear when using several scaled reference systems are 
sketched in Section |6] The system of astronomical units in the 
Newtonian and relativistic frameworks is discussed in Sections 
|2]and[8] respectively. The practical usage of the various scaled 
quantities and also astronomical units in the relativistic context 
is given in Section|9]with the example of extracting the masses 
of the Sun and the Earth from DE405 in SI units. The ques- 
tion of whether the astronomical units of measurements are still 
needed in modern astronomical practice in their current form is 
discussed in SectionfTOl 



2 



S.A. Klioner: Relativistic astronomical time scales 



2. Quantities, their numerical values and units of 
measurements 

In order to discuss the scaling issues let us first clearly distin- 
guish between quantities and their numerical values which ap- 
pear when the quantities are expressed as numbers using some 



units of measurements. According to llSOl (I1993L definition 



1.1), quantity is an attribute of a phenomenon, body or sub- 
stance that may be distinguished qualitatively and determined 
quantitatively. A value (of a quantity) is defined as the mag- 
nitude of a particular quantity generally expressed as a unit 
of measurement multiplied by a number ( llSQ..1993t definition 
1.18). Thus, for any quantity A one has 



A = {A}xx [A]xx, 



(1) 



where {A}xx is the numerical value (a pure number) of quan- 
tity A and {A'lxx is the corresponding unit. Notations [A] and 
[A] for the numerical value and unit of a quantity A, res pec- 
tively , are taken from the international standard ISO 31-0 ( ISOl 
19921). Since in this paper we use several systems of units, the 
subscript gives the name o f the system o f units. Below we use 
index "SI" for the SI units (IBIPM||2006|) and index "A" for the 
system of astronomical units discussed in Section [T] and [8] be- 
low. When a relation is valid with any system of units, as in the 
case of Eq. ([TJ, "XX" is used. 

The of ficial metrol ogical definition of the concept of "unit" 
is given bvliSC 



(119931 



definition 1.7): a unit (of measurement) 
is a particular quantity, defined and adopted by convention, 
with which other quantities of the same kind are compared 
in order to express their magnitudes relative to that quantity. 
Loosely speaking, a unit is a recipe of how an observer can 
realize a specific physical quantity called "unit". The observer 
can then express numerically all other quantities of the same 
kind (those of the same physical dimensionality) by comparing 
them to that specific quantity. The complete set of definitions of 
the concepts o f quantities, values, units, systems of units, etc. 



can be found in llS0l(ll993h . A detailed discussion of these con 
cepts i n the framework of general relativity is given bv lGuinot 
(Il997h . 

Now let us consider two quantities A and B having the same 
physical dimensionality and related by the following formula 
derived in some theoretical way: 



B = FA, 



(2) 



F being a pure number, that is a numerical, dimensionless co- 
efficient. This formula relates quantities A and B irrespective of 
any considerations of units. To get a relation between numeri- 
cal values of A and B one has to use Eq. ^ on both sides of 
Q. In particular, one has 



\B)x 



F{Ah 



(3) 



if and only if one uses the same units for both A and B: [A\xx - 
[B'lxx- Note that one could start this discussion with any kind 
of formula relating A and B with the same conclusion: such a 
formula is also valid for numerical values of the quantities if 
and only if the same units are used for both of them. 



2.1. Units of measurements vs. units of graduation 

Strictly speaking, the concept of "units of measurements" 
can be only applied to measurable (observable) quantities 
(e.g. proper time), but not to non-measurable (i.e. coordinate- 
dependent) qu antities in the framework of general relativity 
(IGuinotil 19971) . For the latter kind of quantities one introduces 
the concept of "units of graduation", which is an alias of "units 
of measurement" for non-measurable quantities. The concept 
of "units of graduation" was introduced to stress that the quan- 
tity under consideration is not measurable so that its "unit" 
cannot be directly realized by some physical measurements. 
However, it seems appropriate to ignore this subtle semantic 
difference in astronomical literature. 

Indeed, let us consider the theoretical formula relating the 
proper time r of an observer with the coordinate time t of some 
relativistic reference system 



T = fit). 



(4) 



Proper time t is a measurable quantity while coordinate time f 
is not. This equation can be derived from the metric tensor of 
the corresponding reference system and from the trajectory of 
the observer in that reference system. Clearly, this is a relation 
between quantities f and r and has nothing to do with units to 
be used at the next step to express these quantities as numbers. 
The same units should be again used for both t and t if rela- 
tion (|4|i is to be valid also for the numerical values of these two 
quantities. If one decides to use the SI second as unit of mea- 
surement for proper time r and if one assumes (|4]i to be valid 
for the numerical values of r and t, the corresponding unit of 
graduation of t is "Sl-second-compatible unit of graduation". 
It is safe, however, to simply call that latter unit of graduation 
"SI second". It seems unnecessary to distinguish "measurable" 
and "non-measurable" quantities at the l evel of units (a discus- 
sion of this point can be found in dGuinot 1997) ): it is sufficient 
to distinguish these two kinds of quantities at the level of their 
physical meaning and properties. In the following we will call 
both units of measurements and units of graduation just "units". 

2.2. On the possibility of scaled units 

Let us note that a linear relation like Eq. ^ could be in prin- 
ciple interpreted as a relation between numerical values of one 
and the same quantity expressed in different units (one quan- 
tity C, two different units related as [C]i = F"' [C]2, so that 
the corresponding numerical values {C)i - F {C}2, and in Eq. 
^ B = {C)i and A = {C)2). However, it is dangerous and con- 
fusing to introduce several units for the same physical dimen- 
sionality (especially, if these units are so close to each other 
that there is a possibility of confusion). The way to introduce 
two different units is against the usual metrological rules (one 
unit for one dimension) and also against the lAU Resolutions 
1991 (Recommendation II) that recommended the use of the 
SI units for all quantities appearing in astronomical coordinate 
systems (in particular, the use of the SI second for all time 
scales). Section |5] below contains a further discussion of the 
topic in connection with the concept of coordinate time scales. 



S.A. Klioner: Relativistic astronomical time scales 



3 



3. Relativistic scaling in the BCRS 



3.2. TDB as a linear function of TCB 



3.1. Dynamical equations in tine BCRS 

Let (f - TCB,x') be the coordinate time and spatial coor- 
dinat es of the Barycentric Celestial Reference Sys tem of the 



lAU (IIAUII200U lRickmanll2001l: ISoffel et al.ll2003h . From the 



BCRS metric tensor one can derive the so-called Einstein- 
Infeld-Hoffmann (EIH) equations of motion of massive bodies 
considered as mass monopoles with masses Ma (capital latin 
indices A, B and C enumerate the bodies): 



Xa 



= -2 



fAB 



j_ y rAB 



V ^ 4 y _^ 3 (rAB ■ Xb) 



CM 1'"'^'^' ^ ''"'^^'^ 

rAB ■ rsc 



Ci'A,B 



\rBc\- 



1 Xa 



■Xb 



BM 



■ 2 ■ Xb - Xa ■ Xa + 4-Xa ■ Xb 



^Xa ■ rAB -^Xb • fAB 



_lZy 



2 ^ \rAB\ 



\rABV 



2j - — + Oifi ), 



\rBC? 



(5) 



and the following equation for the time of light propagation 
between two points x\ and X2 (again for the solar system con- 
sidered as a system of mass monopoles) 



c{t2-tx)^\X2- Xx\ 

2flA 



\riA\ + \r2A\ + k2i 
\r2A\ + \riA\ - k2i 



(6) 



where ^a - G Ma, G is the Newtonian gravitational constant, c 
is the light velocity, Xa is the position of body A, a dot denotes 
the time derivative with respect to TCB, and for any two indices 
fAB = Xa- Xb- Quantities ^a - G Ma are called "mass param- 
eters" below in order to clearly distinguish them from masses 
Ma and the gravitational constant G. Eq. (|5]l was derived as 
early as in 1917 by Lorentz and Droste and then re-derived 
by Einstein, Infeld and Hoff mann from a more general point 
of view (see Section VII.C of lDamour. Soffel & Xul(IT991,) for 
a detailed history of these equations). The EIH equations have 
been used to construct accurate solar system ephemerides start- 
ing from the middle of the 1970s. Eq. ^ describes the rel- 
ativistic Shapiro time delay which is also well known since 
the 1960s and is widely used for astronomical data modeling. 
Note, that up to some theoretical improvements, the BCRS was 
known akeady in the 1930s and even earlier. The lAU has only 
officially fixed the status quo in the Resolution B1.3 (2000). 
In both equations above the coordinate time t - TCB of the 
BCRS is used. It is TCB (and not TDB or any other time scale) 
that was used since 1917 in all theoretical works involving the 
aforementioned equations and the underlying relativistic refer- 
ence system. 



For the reasons of practical convenience one often uses the so- 
called f — TDB which is a linear function of f = TCB: 



f - Ft + t\ 



0' 



(7) 



with F — \ - Lb and are defining constants adopted by the 
lAul (l2006h . The constants are fixed here in such a way that 
TDB, evaluated at the geocenter, remains as close as possible to 
TT (see Section|4]below). In particular, the mean rate of TDB 
coincides with the mean rate of TT when the transformation 
between TDB and TT is evaluated at the geocenter. The mean 
rate of TT in its turn is equal, to a high level of accuracy, to 
the mean rate of the proper time of an observer situated on the 
rotating geoid. TT is directly available to the Earth-bound ob- 
servers through TAT UTC or any other realizations of TT. The 
difference between TDB and TT does not exceed 0.002 s and 
can be neglected for many applications. These circumstances 
and lower risk of a damage if TDB is confused with TT (com- 
pared to possible damages of confusing TCB and TT with their 
linear drift of about 0.5 s per year) are the arguments usually 
put forward in favor of TDB. Here we use the new definition 
of TDB recently adopted by the llAUl (l2006l) . The original lAU 
wording given in 1976 defines TDB so "that there be only pe- 
riodic variations between these time-scales" (TT and TDB). 
This definition is known to be fundamentally flawed since the 
resulting TDB is not a linear function of TC B and cannot be 
used with usual d ynamical equations (|5]l-(|6]l ( StandishI 1998 : 
Soffel et al.lE003h . 



Two realizations of the TDB widely used in practice are 
given by the analyti cal formulas relating TDB and TT given by 
Mover (fl981h and iFairhead & BretagnonI (Il990l) . The former 
has a lower accuracy of about 20 fis and contains only periodic 
terms as an attempt to adhere to the lA U (1976) description of 
TDB. The more accurate formulation of lFairhead & Bretagnon 
(1990) contains many non -periodic (polynomial and mixed) 
terms. This demonstrates that retaining only periodic terms in 
the transformation between TDB and TT is only possible as a 
numerical approximation for lower accuracies and shorter time 
spans for which the analytical formulas should be valid. 

Another time scale, very similar to TDB and also lin- 
early related to TCB, was described by IStandishI ( 119981) and 
is called Tepi,. The subtle difference between T^ph and TDB 
lies in the way the constants in (O are chosen. For TDB 
the constant Lb = 1.550519768 x 10"^ in (EJ is a defining 
one while for T^pi, the constant F is different for different 
ephemerides and implicitly defined by the transformation be- 
tween TT and Tgpi, used during the construction of each partic- 
ular eph emeris. The adopted rPZj value of Lb is based on the 
work of Irwin & Fukushima (Il999h and iHarada & Fukushimal 
(l2003l) . and on the lAU Resolution B1.9 (2000) defining TT. 
The additive constant in (|7]i plays no role for the purposes of 
this paper and will not be discussed here. 
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3.3. Scaling of spatial coordinates and mass 
parameters 

If one uses f* = TDB instead of TCB it is natural also to intro- 
duce scaled spatial coordinates x* and scaled mass parameters 
fi* for each body as 



X* - Fx, 



(8) 
(9) 



These additional scalings allow one to retain exactly the 
same form of the principal dynamical equations (|5]l and (|6]l. 
Quantities x and fi are called TCB-compatible quantities (or 
simply TCB quantities) representing spatial coordinates and 
mass parameters. Quantities jc* and //* can be called TDB- 
compatible quantities (or simply TDB quantities). Let us make 
here several comments: 

(1) Physical mass of a body corresponds to ju (not to yU*). 
Quantity fi does not depend on the kind of experiments used 
to get it and on where the observer measuring it is situ- 
ated. Moreover, fi represents also the mass parameter of the 
corresponding body in the Ge ocent ric Celestial Reference 
Syste m (OCRS) of the lAU dlAUl EoOl: Rickman 2001; 
ISoffel etalJ 120031) . On the other hand, the scaling factor 
between /u and fj* is related to the fact that most of accu- 
rate observations were until now performed from the sur- 
face of [rotating] Earth. This made TT (TAI or UTC) con- 
venient or even natural to parametrize observations. This 
will change as soon as sufficiently large amount of high- 
accuracy observations is performed from space vehicles. 
Therefore, //* can only be considered as some ad hoc pa- 
rameter convenient from some practical point of view. 

(2) It is confusing to believe that t, f, x, x*, fi, fj* are not 
quantities, but just values in different units: f, x and /u are 
numerical values expressed "in SI units", and f*, x* and 
fj.* are "in TDB units". As discussed in Section l272l Eqs. 
(Q, (HI) and (|9]l are relations between six distinct quanti- 
ties, and the question of units has not been discussed at all 
here. One can, for example, consider (|7]l and ([8]) as rela- 
tivistic coordinate transformations from {t,x) to {t*,x*) in- 
troducing another reference system BCRS* distinct from 
BCRS. On the other hand, some "TDB units", distinct to 
the SI units, imply, in particular, that the "TDB second" is 
no longer SI second (and the number "9192631770 " as ap- 
pears in the definition of SI second (BIPM 2006f) should 
be explicitly changed to some other number in the defini- 
tion of that "TDB second"). Those non-SI units have never 
been defined or discussed seriously. Moreover, since there 
exists a consensus, also enforced by the lAU Resolutions 
(1991), that TT, TCG, TCB and related quantities are all 
measured in the SI units (or in "Sl-compatible units of grad- 
uation"; see Section ITTI )). it is quite natural to extent this 
rule to TDB and TDB-compatible quantities and to make 
the semantics symmetric and clear In Section |5]below one 
can find additional arguments against the concept of "TDB 
units". 



3.4. Further implications of the three scalings 

Eq. ^ relating x and x* is valid for any distance used simulta- 
neously with TCB and TDB. In particular, the TCB-compatible 
semi-major axis a of a planet is related to the corresponding 
TDB-compatible semi-major axis a* as 



a* — Fa. 



(10) 



In the same way Eq. ^ also implies corresponding relations 
between time intervals. In particular, the TCB-compatible and 
TDB-compatible orbital periods of a planet are related as 



FP, 



(11) 



and the corresponding mean motions (n = 2nlP, n* - 2nlP*) 
as 



F-^n. 



(12) 



CoiTespondingly, the third Keplerian law for a massless particle 
moving in the field of a central body reads 



3 2 

a n — ji 



for the rCB-compatible a, n and yu, and 



*3 *2 » 

a n - ju 



(13) 



(14) 



for the TDB-compatible ones. 

As discussed in Section |2] relations (l7Ti-(fl4]i are also valid 
for numerical values of the corresponding quantities if the same 
units are used for the quantities appearing on both sides of these 
equations. If the units used for TCB- and TDB-compatible 
quantities are the same one has 



{tlxx ^ F lt]xx, 
{x*]xx ^F{x}xx, 
il^lxx ^FWxX' 
{a*}xx = F{a}xx, 
{nxx^F{P}xx, 
{n}xx = F'^ {n}xx, 
ffllix {nfxx = Mxx' 
{afxxinfxx = 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



where subscript "XX" denotes the name of any chosen system 
of units. Those "same" units could be the SI units (SI seconds 
and SI meters) as recommended by the lAU, but also any other 
system of units: for example, astronomical units that are widely 
used in astronomy during the last two centuries. 

4. Relativistic scaling in the GCRS 

Let us consider the Geocentric Celestial Reference System 
(GCRS) of the lAU with coordinates (T = TCG, X). For the 
reasons discussed above for TDB, it is often convenient to in- 
troduce a scaled version of TCG called T** - TT. For current 
clock accuracies the mean rate of TT is indistinguishable from 
to the mean rate of the proper time of an observer situated at 
the rotating geoid. The difference comes from the tidal effects 
and does not exceed 10 in the rate and 1 ps in the amplitude 
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of periodic terms. Again the scaling of time coordinate makes 
it convenient to introduce the scaled versions of spatial coordi- 
nates and mass parameters: 

T" ^LT, (23) 
Z"=LZ, (24) 

(25) 

with L = 1 - Lr„ L^, = 6.969290134 x 10"'" being a defining 
constant (lAu" 2001). As discussed above ^ is the same in both 
GCRS and BCRS. Here again we argue that one should speak 
about six independent quantities (three TCG-compatible quan- 
tities T, X and ji and three TT-compatible ones T**, X** and 
//**) without any reference to units. The expression "TT units" 
is confusing for the same reasons as "TDB units" and should 
be avoided. Numerical values of these six quantities expressed 
in the same units are scaled in the same way as the quantities 
themselves. For example, in the SI units one has 





(26) 


[X"}si^L{X]si, 


(27) 




(28) 



5. On the concept of coordinate time scales 

The question if and in which sense the SI second can be used 
with the coordinate time scales is closely related to the under- 
standing of what coordinate time scales are. Although the con- 
cept of a coordinate time is crystal clear for people trained in 
relativity, coordinate time scales maybe sometimes very con- 
fusing for people using "Newtonian common sense". In the lit- 
erature one can sometimes meet wrong statements about rela- 
tivistic coordinate time scales. Among these wrong statements 
one can find: (1) TCB is the time in the barycenter of the so- 
lar system, (2) TCG is the time at the geocenter, (3) TT is the 
time on the rotating geoid, (4) an ideal clock put in these three 
locations would keep TCB, TCG and TT, respectively, (5) for 
TDB no location could be found where an ideal clock would 
keep it and this implies "TDB seconds". All these statements 
originate in an incomplete or inconsistent knowledge of rela- 
tivity resulting in a latent yearning to save Newtonian absolute 
time and at least some of its nice features. Let us try to depict 
the role of coordinate time scales. 

Let us consider first proper time of an observer. Since the 
SI second pla ys an import ant role for this discussion let us cite 
its definition (lBIPMll2006h : 

The second is the duration of 9 192 631 770 periods of the 
radiation corresponding to the transition between the two hy- 
perfine levels of the ground state of the caesium 133 atom. 

It is important to understand that, in full agreement with 
general relativity, this definition contains no indication of any 
specific location, gravitational potential or state of motion of 
the observer realizing the SI second. This means that this def- 
inition is a recipe how any observer can realize the SI second. 
For all observers this recipe is the same and in this sense the 
SI second is the same for all observers. This means also that an 
observer has no chance to notice his motion and position look- 
ing only at the readings of its clock. Therefore, the SI second 
by itself is the unit of proper time. 



The differences between proper times of two different ob- 
servers can only be noticed when a comparison procedure for 
two clocks having different trajectories is established. If these 
observers are located at the same place at the moments of com- 
parison (as, e.g., in the twin paradox), the comparison proce- 
dure is obvious. Otherwise (which is the typical case) it in- 
volves some (arbitrary) relativistic reference system which de- 
fines the coordinate simultaneity. The concept of simultaneity 
is an indispensable part of any clock comparison algorithm. 
The coordinate simultaneity is the only logical possi bility to 
save t he concept of simu ltaneity in general relativity jKlionei 
19921: IPetit& Wold l2005h two events are called simultaneous 



if and only if the chosen coordinate time has the same value for 
both of them. 

Let us turn to the coordinate time scales TCB, TDB, TCG 
and TT. These time scales are part of the mathematical model 
of space-time used in general relativity. The mathematical 
model of spaced-time is called reference system and represent 
a 4-dimensional chart allowing one to assign four numerical 
labels for each space-time event. Three of these four labels are 
called spatial coordinates and the fourth label is called coordi- 
nate time. Coordinate time scales TCB, TDB, TCG and TT are 
defined for any space-time event within solar system and far 
beyond. All these coordinate time scales are coordinates and, 
therefore, cannot be directly measured. They can only be com- 
puted from the readings of some real clock(s). For this compu- 
tation one should use a theoretical relation between the proper 
and coordinate time scales that follows from the basic princi- 
ples of general relativity. That theoretical relation involves cer- 
tain model of the solar system: the trajectory of the observer, 
the trajectories of the massive bodies, their mass parameters, 
etc. For TCB and TDB this model is given in the BCRS, and 
for TCG and TT in the GCRS. 

Sometimes, especially for didactic reasons, it is useful to 
consider a special imaginary observer (that is, a special trajec- 
tory for an imaginary observer), the proper time of which coin- 
cides with the considered coordinate time along the observer's 
trajectory. For example, for TCB such an observer is situated 
infinitely far from the solar system (so that the gravitational po- 
tential of the solar system vanishes at his location) and is at rest 
relative to the solar system barycenter Certainly one can find 
also an analogous observer for TDB: take the same observer 
as for TCB, but moving with a constant velocity so that the 
Lorentz time dilation exactly compensates the rate difference 
between TCB and TDB. Note however, that such observers are 
only useful as an illustration and their existence should be over- 
estimated. First, these observer does not help to define TCB or 
TDB: the proper time of an observer is defined only on his tra- 
jectory, while TCB and TDB are both defined everywhere in 
the solar system. Second, these observers do not help to relate 
a real clock moving within the solar system to TCB and TDB. 
Third, for TCG and TT it is not possible to find such imaginary 
observers. One often argues that TT "is defined on the rotating 
geoid". This is not true since TT is a coordinate time scale and 
is defined for any event in the solar system. However, the rela- 
tion between the proper time of an observer and TT shows that 
the proper time of an observer situated on the rotating geoid 
is close to TT computed along his trajectory, but only up to 
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terms of order 10 and up to periodic terms of an amplitude 
of about 1 ps. It is just a close agreement for a particular tra- 
jectory and no more than that. TT can be computed along any 
other trajectory and can be related to the proper time of any 
other observer in the solar system. Note also that the constant 
Lc in the definition of TT is decoupled from the geoid and will 
not be changed when the definition of the geoid is improved. 

From the theoretical point of view the situation with the 
pair TCG-TT is completely symmetric with the situation with 
the pair TCB-TDB: TT is a scaled version of TCG, while TDB 
is a scaled version of TCB. Both scalings have no physical 
meaning, but were chosen from the considerations of conve- 
nience: to make the difference between the proper time of an 
observer on the rotating geoid and these two coordinate time 
scales evaluated along his trajectory as small as possible. 

From the metrological point of view, the theoretical rela- 
tions between the proper time and the coordinate times are re- 
lations between quantities and are independent of the choice of 
units. The same is true for the theoretical relations between the 
coordinate time scales themselves. This means that any units 
can be used to convert these quantities into numerical values. 
The SI second is usually used for the proper time of any ob- 
server. If the theoretical formulas related the proper time and 
coordinate times are used to compute the values of the cor- 
responding coordinate time starting from the values of proper 
time in SI seconds, the values of the coordinate time are ex- 
pressed also in SI seconds (see. Section ITTT i. The same is true 
for the theoretical formulas relating coordinate time scales with 
each other. In this situation some non-SI "TDB units" are com- 
pletely artificial and unnecessary. 

6. Mixing scaled BCRS and scaled GCRS 

The scaling of BCRS and GCRS is obvious and simple to man- 
age if only one of these reference systems is used. In practice, 
however, relativistic models often involve quantities defined in 
both reference systems. Good examples here are models for 
VLBI and LLR, but it is also the case for virtually all kinds of 
observations. For example, VLBI model contains station coor- 
dinates and Earth orientation parameters defined in the GCRS 
while the positions of sources and solar system bodies (e.g. the 
Earth and the Sun) are defined in the BCRS. These "mixed" 
models are not invariant under the scalings (l7]i-(|9]l and (l23Tl- 
$25[ . As a result the coefficients Lb and Lq (and the constant 
Lc = 1 - (1 - Lb) /{I - Lc) ~ Lb- Lg) expUcitly appear in the 



standard VLBI model (IIERS Conventions 1120031 Chapter 11) 



This makes the models less transparent conceptually and more 
difficult to understand and maintain. It should be stressed that 
the scalings (and the corresponding coefficients) represent non- 
physical, conventional changes of the BCRS and the GCRS and 
do not appear in normal relativistic considerations. 

Moreover, with the increasing importance of spacecraft ob- 
servations, the number of coordinate systems we have to deal 
with has proliferated. For example, in order to study local 
physics (e.g. rotational motion) of Mars, Moon or Mercury one 
could introduce GCRS-like planetocentric reference systems in 
the vicinity of each planet. In particular, these local planeto- 
centric reference systems introduce their own coordinate time 



scales. Consistent adherence to the idea of scaled coordinate 
times having the same rate as TT at the centers of mass of 
these planets would require special scaling factors for each of 
these reference systems. This would make the data reduction 
schemes disastrously complicated and obscure. 

In principle, it would be cleaner from the point of view of 
theoretical purity and consistency not to introduce these scal- 
ings at all. It is however clear that the considerations of conve- 
nience and a kind of tradition weigh against the full use of the 
original non-scaled versions of the BCRS and GCRS. 

7. The system of astronomical units in Newtonian 
framework 

The reason to introduce astronomical units of measurements in 
the 19th century was the fact that the accuracy of positional 
(angular) observations was much higher than the accuracy of 
determination of distances (e.g. solar parallax). Before the in- 
vention of radar and laser ranging and related techniques it was 
much easier to measure the period of motion of a planet than to 
determine the distance to that planet from the Sun or from the 
Earth (only a kind of geometrical triangulation could be used: 
e.g. observations of Venus transits or of Eros in its close ap- 
proach to the Earth). For that reason, solar system ephemerides 
have always been first constructed in the so-called astronom- 
ical units to use the full precision of positional observations 
and only later (and only if necessary) were they converted into 
other units directly available in a laboratory (e.g., metric units). 
The precision of that last conversion could be [much] lower 
than the precision of the ephemeris in astronomical units. The 
ephemeris in astronomical units is sufficient, however, to pre- 
dict angular positions of the bodies on the sky. 

Let us first forget about relativity and consider the classi- 
cal Newtonian situation. The system of astronomical units con- 
sists of three units: one for time f, one for mass M and one for 
length X. From now on we designate these astronomical units 
as [f]^, {M}a ™d [jc]^, while the corresponding SI units are 
[f]^/ = second, [Al]^/ = kilogram and {x\si = meter. The cor- 
responding numerical values in astronomical units are denoted 
as {t]^, {M]a and {x)^, and in SI units {f)s/, {Al}^/ and {x]^!- 

The astronomical unit of time is the day. The day is directly 
related to the SI second: 



[t]A = day [t]si. 



(29) 



where d - 86400 is a pure number The astronomical unit of 
mass is fixed to coincide with the "solar mass" (SM) 



[M]a = SM = q' [M] 



(30) 



where a is a pure number giving the solar mass in SI kilo- 
grams. The value of a should be determined from observa- 
tions (see below). The astronomical unit of length [x]^ is called 
"Astronomical Unit" (AU) 



[x]a ^AU^x [x]si. 



(31) 



where x is the number of SI meters in one AU. The AU is de- 
fined in a tricky way with no relation to x- First, one fixes the 
value of the Newtonian gravitational constant G expressed in 
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astronomical units to coincide with the value determined by 
Gauss in 1809 from a series of observations available to him. 
For historical reasons that value is used up to now as a defin- 
ing constant in the definition of the system of astronomical 
units (k is the well-known and widely-used standard notation 
for Vicjl): 



{GU =k^ = 0.01720209895^ 



0.000295912208285591 1025. 



(32) 



Clearly, in any system of units the dimensionality of G is 
[jcJxx^ [^Ixx ' (^nd in particular, in astronomical units 
the dimensionality is AU'^day^^SM '). The AU is then defined 
to be the unit of length with which the gravitational constant 
G takes the numerical value ( |32] |. This definition of the as- 
tronomical unit was adopted by the lAU in 1938. One can 
also say that the AU is the semi-major axis of the [hypothet- 
ical] orbit of a massless particle which has exactly a period of 
^ * 365.256898326328 . . . days (astronomical units of time) 
in the framework of unperturbed Keplerian motion around the 



Sun having the mass of 1 SM dBrouwer & Clemencd 11961 



Standishll2005al) . Kepler's third law gives 



{a]l {n}l = {GU {MU, 



(33) 



where {a]^ and {n]^^ are the numerical values of the semi-major 
axis and mean motion of a Keplerian orbit expressed in astro- 
nomical units of length and time, respectively, and {M}^ is the 
mass of the central body in the astronomical units of mass SM. 
Normally, in the classical Newtonian case for the mass of the 
Sun Mq one can just put {AIq}^ = 1. Note, however, that SM 
is a unit, while Aie is a quantity. Moreover, SM must not co- 
incide with the real physical mass Me of the Sun, especially 
since the latter is time-dependent (see SectionfTOli. 

The system of astronomical units is defined by four num- 
bers d , a, x k. In modern astronomical practice CStandish 
2005 at) the value of;^' is determined from the whole set of avail- 
able observations (various ranging observations that measure 
distances directly in SI units play here a crucial role). Then, 
comparing ( |33] ) and ( 1211 1. the numerical value of fi - G Mq in 
SI units can be computed as 



SI 



k^{MeUx'd-\ 



(34) 



The mass of the Sun in kg can then be computed by using the 
SI value for G ({G}si = 6.674 . . . x 10"" m^ s^^ kg"'), but this 
last step is not important for precise work. 

Using the relations between the astronomical and SI units 
one can write the following relations between numerical values 
of time f, distances (positions) x and mass parameters fi 

{tU = d-' {t}s„ (35) 

{xU^X'Uxhi, (36) 

{fiU = d^X'' {f^}sn (37) 

and for the period P, mean motion n and semi-major axis a of 
an orbit 



{P}a 
Ma 



: {Pi 



Sh 



d [n 



SI- 



(38) 
(39) 
(40) 



8. The system of astronomical units in tlie 
relativistic framework 

Up to recently, only TDB was used as independent time ar- 
gument of modern ephemeris. In connection with efforts to 
construct new ephemerides with TCB (or to re-parametrize 
old ones) the system of astronomical units in the relativistic 
framework has been considered recently by several authors 
(fStandis h .19951: Urumberg & Simon 2004; Standish 2 005bl: 
|Pitievall2005bl) . Let us interpret here all the formulas in the 
previous Section as formulas relating TCB-compatible quan- 
tities. As we discussed in Section[3]the TDB-compatible quan- 
tities are related to the corresponding TCB-compatible ones by 
a relativistic scaling. Then one can introduce another "TDB- 
compatible" system of astronomical units (designated as "A*"): 



[t]A* = day* = d* [tjsi, 
[Mh. = SM* = a* [Mhi, 
[x]a, = AU* = X* Msi, 



(41) 
(42) 
(43) 
(44) 



Let us first consider the four numbers d*, a*, x* and k* as ar- 
bitrary (totally independent of the corresponding four numbers 
d, a, X and k defining the TCB-compatible system of astro- 
nomical units) and write down the relations between numer- 
ical values of TCB-compatible quantities expressed in TCB- 
compatible astronomical units (e.g. (ju)^) and TDB-compatible 
quantities expressed in TDB-compatible astronomical units 
(e.g. {jU*UJ: 



(45) 
(46) 
(47) 



itlA.-F\-\ {tU, 



F ^- 



X 



These relations should be compared to the corresponding 
relations ([TsTl-lfTTb in SI units. Considering that {;U*}^» - 
{k*f {AIq}^^ and \pi}^ - k^ {Mq}a one can rewrite Eq. WT\ 

as 



kHM^}A 



(48) 



This is the only constraint on the involved constants. Starting 
from this relation one can suggest many different ways to define 
both TCB- and rZJB-compatible systems of astronomical con- 
stants. One reasonable additional consideration is that a "day" 
is defined to be 86400 seconds in any time scale (in TCB, TDB, 
TCG, TT or proper time of an observer) as was recently agreed 
by the lAU Working Group on Nomenclature in Fundamental 
Astronomy (Capital n3l2005 ). This means that the physical du- 
ration of a day depends on the time scale used. Therefore, "day" 
is defined by the conversion factor 86400 and it is natural to put 
d* - d - 86400. Considering this one has at least two choices: 

I. One can require that the solar mass is equal to 1 in corre- 
sponding astronomical units in both cases and that the con- 
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Slants k* and k are equal as well (IStandishll995h . This gives 
{a**}/!* - Wa - ^iid together with d* - d leads to 



This, in turn, gives 

{/^*!a. = K^Aa- 



(49) 



(50) 
(51) 
(52) 



9. Numerical example: mass parameters, 
coordinates and velocities from DE405 

Let us illustrate how to extract numerical values of planetary 
mass parameters in SI units from the existing ephemerides con- 
structed using TDB and the corresponding system of astro- 
nomical units, using as a specific example JPL's DE405. In the 
header of DE405 one finds the following TDB-compatible nu- 
merical values: 



Note the unusual scaling laws of distances expressed in 
astronomical units and the astronomical unit itself in this 
case. These scaling laws have produced akeady a lot of 

confusion in the liter ature. 

II. Another possibility ( Brumberg & Simon 2004; StandishI 
2005bl: lPitjevall2005 b) is to retain the scaling laws of time. 



= 1.49597870691x 10", 

= 2.95912208285591 1025 X 10" 



(57) 
(58) 



distance and mass in astronomical units (|45]|-(|47]| in ex- 
actly the same form as in SI units (fT5ll-(fT7]l and put 



The latter number is just the Gaussian value of ]^ quoted in 
( |32] |. This allows us to find from ([34b 

= {j"g1a. Or*)^ 86400^2 = 1.32712440018 x 10^°, (59) 
and for the TCB-compatible mass parameter 



{//ols^ = [^1%]^^ = 1.32712442076 X 10 



20 



(60) 



which together with d* - d gives 

{f*)A. ^^'{flA, 

{x*\a* ^F{x)f,, 
{/!a. ^F\ii\,^ . 



(53) 



(54) 
(55) 
(56) 



That latt er value for the sol a r mass parameter can be found, 
e.g., in llERS Conventions! (l2003b . The TDB- and TCB- 



compatible masses of planets can be found from their TDB- 
compatible masses in astronomical units given in the header of 
DE405 in the same way. For example, for the Earth one has 



{//^l^^ = 3.98600432889 X 10 



14 



The only "unusual" consequence of this choice is that 
(A;*)^|A40|^ it k^ [Mq]j^. Since modern ephemerides 



■els I ■ 



F-' ll^l 



3.98600439069 X 10' 



The corresponding TT-compatible value is then 



constructed with rPfi-compatib le quantities use {^;*}^^ = ^{^(,3)5, = 3.98600438792 x 10 



(r)2(M*}^^ = 0.01720209895^ (IStandishI Eo05bl) this 
means that for TCB-compatible units one has k^ {AIqIa - 
F-i 0.01720209895^ ~ 0.00029591221287376846 . . .. 

The second choice seems to be more natural. Let us note 
finally, that contrarily to what can be inferred from some pub- 
lications, the definitions of astronomical units by no means 
influence the relations between numerical TCB- and TDB- 
compatible quantities (e.g. yu and yu*) in SI units: they remain 
to be defined as shown in Eqs. (fT5])-(l20l). 

It is unclear what role these "scaled" relativistic astronom- 
ical units could play for new solar system ephemerides: astro- 
nomical units are just units and one can use any definitions of 
them as long as the definitions are known. It makes no sense 
just to reformulate the same process to produce ephemerides 
with TCB (instead of TDB) and TCB-compatible astronomi- 
cal units (instead of TDB-compatible astronomical units): the 
results after corresponding re-scaling must be identical to the 
TDB-compatible ones, provided that all the scalings are per- 
formed in a consistent way. If all scaling factors appearing in 
the process of ephemeris development were applied correctly 
one can claim that with the same level of confidence one can 
just re-scale an existing ephemeris constructed in TDB into 
TCB according to equations given above. The question of con- 
sistent use of relativistic time scales in the process of construct- 
ing new so lar system eph emerides will be considered in detail 
elsewhere ( Klionei|[2007 ). 



(61) 
(62) 

(63) 



Let us also note that the mass parameters of the major plan- 
ets, except for Pluto, are all based currently on observations of 
spacecraft motions in the vicinity of the corresponding planet. 
This means that an additional re-scaling procedure should be 
performed between the mass parameters used for different ap- 
plications (e.g., for the Earth {/^e j^^ used for SLR, {jt'e}^/ f^^" 
the planetary ephemerides etc.). 

Let us now turn to positions and velocities. The DE data 
in the distribution gives the numerical values of the TDB- 
compatible spatial coordinates x* in SI units, that is {x*}s[ 
parametrized by t* (precisely speaking the coordinates are 
given in kilometers, not in meters, but it plays here no role). 
If the TCB -compatible positions x are desired in SI units they 
can be computed as (cf. Eq. ( fTSI l) 



isi 



isi- 



(64) 



The TCB-compatible velocity coincides with the TDB- 
compatible one v* = dx*/dt* = dx/dt = v. The standard 
JPL software can be asked to output either directly {x*]si or 
the TDB-compatible position JT* in TDB-compatible astronom- 
ical units, that is {jr*)^,. The latter value is computed by divid- 
ing {x*]si by X* from the header of the ephemerides (see, Eq. 
(I57I1). A similar procedure applies to the velocity v*. Note that 
it is only;^^*, and therefore, only TDB-compatible astronomical 
units that can be considered as a part of the DE ephemerides. 
The choice of the TCB-compatible astronomical units (and in 
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particular the value ofx) is by no means influenced by the DE 
ephemerides or by the procedures used during their develop- 
ment. Let us stress again that this choice has no influence on 
the relation between the values in SI units. 



10. Do we need astronomical units in their current 
form? 

It is not clear if astronomical units should be further used to 
construct future ephemerides. The main reason for astronom- 
ical units - much higher accuracy of angular (positional) ob- 
servations compared to distance measurements - does not exist 
any longer Considering the subtleties with astronomical units 
in the relativistic framework one can find it more advantageous 
either to avoid astronomical units at all or for reasons of histor- 
ical continuity to convert them into defined units by fixing x as 
it was done with the day (d = 86400) and with the SI second. 

One more argument against the system of astronomical 
units in its current form is that the physical mass of the 
Sun is, in principle, not constant, but d ecreasing at a rate 
o f ~ 10~" solar masses per century (INoerdUngerl 1 1997 ; 



Krasinski & Brumbe rg 2004) just because of the Solar radia- 
tion. Up to now the dynamical consequences of this change 
were below the accuracy of observations, but one can expect 
that in the near future astronomical measurements in the inner 
solar system will reach a level of accuracy where the effects 
of changing solar mass (secular acceleration in the mean lon- 
gitudes of the planets) will become observable. For example, 
Pitieval (12005 ah gives the accuracy of the determination of GIG 
as 5x10 per century. This is the precision of the claim that no 
secular accelerations in the mean longitudes of the inner plan- 
ets are observable. On the other hand, a linear change of the 
mass of the Sun has the same consequences for astronomical 
observations as a linear change of G. Thus, in the near future 
we will have to decide if we want to live with time-dependent 
units of length, fix some epoch to define the Astronomical Unit, 
avoid astronomical units in precise work, or, preferably, make 
the Astronomical Unit to be a defined constant by fixing y for 
historical continuity. 

If the AU is fixed in SI meters, the mass of the Sun or, 
more precisely, fi^ - GMq should be fitted from observations 
together with masses of other planets, while the AU plays the 
same role of a "convenient" unit as kilometer or mile. It seems 
to be even more reasonable since in modern practice the masses 
of the planets are often determined by other kinds of observa- 
tions that deliver yu directly in SI units. For example , the curren t 
best value for fig, is delivered by SLR dOrotenl 1 999t lRiedl2005h 
with no relation to astronomical units. 
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